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obtained explicitely. 
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In recent years much attention has been focused on an issue of the criticality 
and phase transitions in ionic fluids. For reviews of the experimental and theoretical 
situation, see Ref. PQ- |3j- One of the most widely used theoretical models of Coulomb 
systems is the so-called "restricted primitive model" (RPM), in which the ionic fluid is 
modelled as an electroneutral binary mixture of charged hard spheres of equal diameter 
immersed in a structureless dielectric continuum. Early studies [4 established that the 
model has a gas-liquid critical point (CP), the location of which was determined by 
recent simulations Due to the controversial experimental findings, the critical 

behavior of the RPM has also been under active debates and strong evidence for an 
Ising universal class has been found by computer simulation [3 IE] as well as by a recent 
theoretical study jU]. More recently, the effects of size- and/or charge-asymmetry on the 
phase diagram have been studied theoretically (THj- [12] an d by simulation [15^ fl6| . 

The more relevant issue concerns the behaviour of the two-body charge-charge 
correlation function God{i^) near the CP where the density fluctuations diverge strongly. 
Recently Aqua and Fisher [TH] using a class of exactly soluble spherical models studied 
both symmetric and nonsymmetric version of 1 : 1 ionic models. They showed that 
in the former case the two-point charge correlations remain of short range and obey 
the Stillinger-Lovett (SL) rule near and at the critical point. Otherwise they found the 
divergence of the charge correlation length to be precisely the same as the divergence of 
the density correlation length. They also found that the SL rule fails at criticality. 

In this letter we address an issue of the behavior of the correlation functions of 
an ionic fluid far and near the CP. For this purpose we consider a charge- and size- 
asymmetric continuous model of an ionic fluid. According to [T3| the Fourier transform 
of the charge-charge correlation function of the model Gqq(/c) can be presented as 



where the first vanishing term results from electroneutrality (zeroth-moment SL 
condition) while the second term yields the second-moment SL rule with £z,i — £d 
(£d is the Debye length). We intend to answer the question whether the SL rules are 
satisfied in the cases of both non-symmetrical and symmetrical ionic fluids. 

In order to fix our notations we first consider the case of a fluid made of two species 
a = 1, 2. In the grand-canonical (GC) ensemble an equilibrium state is characterized by 
the inverse temperature (3 = (ks and T are Boltzmann constant and temperature, 

respectively) and two dimensionless local chemical potentials z/ a (r) = (3{[i a — tp a (r)), 
where /i Q is the ordinary chemical potential and ip a { v ) is some external field acting on 
particles. Recall that the GC partition function 5 [i>i, z/ 2 ] is a log-convex functional of 
the v a . The truncated two-body correlation functions are defined as 




(1) 



p>2 



G a p [ui,v*] (1,2) 



5 2 In 5 [i/i,i/ 2 ] 
8v a {l)8v p {2) ■ 



(2) 



AA[pi,p 2 ] = sup <^ ^2 (PaWa) -lnS [vi,v 2 ] 

V1 '^ [ a 
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The Legendre transform of In H defined as 

(3) 

gives us the free energy. In eq. (J3j) Dirac notations have been used, i.e. 

(paWa) = / d 3 r p a (r)u a (r) , (4) 

where f2 denotes the volume of the system. PA [pi, p 2 ) is a convex functional of the local 
densities p a (l) and it is the Legendre transform of lnS [v\, z/ 2 ], i.e. 

lnS [u u v 2 \ = sup (PaWa) ~ PA [pi,p 2 ] \ ■ (5) 

pi>/» { a J 

/3A[pi,p2] is the generator of direct correlation functions. In particular the two-body 
direct correlation functions are defined as 

g^k^(l,2) = - WM2) - (6) 

The functions G and C are related by the Ornstein-Zernike (OZ) equation 

C Q/3 (1, 2)6^(2, 3) = -5^(1, 3) , (7) 

where summation and integration over the repeated mute indices are meant. 

We now consider the case of ionic mixtures, particles of species a = 1 are supposed 
to bear a charge +q and those of species a = 2 bear an opposite charge — Zq (Z > 0). 
It is obviously convenient to introduce the total number and charge densities, pn and 
Pq, by the relations 







AT> '» (8) 




with 



For a homogeneous fluid we certainly have the charge neutrality condition pg = (pg) = 
where the hat denotes the microscopic charge density and the brackets (...) denote a 
GC average. The (truncated) correlations of pn and pq will be defined as Gab(1,2) = 
(p A {l)p B {2)) - (p A (l)} (p B (2)>, where A(B) = N,Q. 

It follows from equations given above that we can introduce two symmetric matrices 
of correlation functions, namely 

G = \ < ^ 11 ] G = ( ^ NN ^ NQ ) (10) 

I G 2 i G 22 J \ G QN Gqq J 

which are related by the matricial equations 

G = M- 1 GM- 1 (11) 
G = MGM. (12) 
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(13) 
(14) 

(15) 



Similarly we introduce the C and the C as 

C = M CM 
C = M~ 1 CM' 1 , 

in such a way that the form of OZ equation is preserved, i.e. 

C(1,3)G(3,2) = -U 5(1,2), 

where U is the unit matrix. 

We are interested in the behavior of Gab{^) in the critical region at large r or 
equivalently at small k in Fourier space. This study will be made in the framework of 
the random phase approximation (RPA) defined by the closure relation 

C = C HS - P$ - /3$ c • (16) 

In eq. f!16|) C*hs is the matrix of the direct correlation functions of some reference system 
chosen here to be a hard sphere fluid (HS) mixture (with diameters a a , a = 1,2). The 
matrix $ is built from short range pair interactions and, finally, $ c denotes the matrix 
of Coulomb interactions. With these notations, the OZ equation ()15j) can be rewritten 
in Fourier space as 

G NN (k) G NQ (k) 
G Q N(k) G QQ (k) 

C* s N (k) - (3^ NN (k) C™(k) - (3^ NQ (k) 



(17) 



C™(k) - 0* NQ (k) Cgg(fc) - P® QQ (k) + 



47T/V 

k 2 



where 



®QQ(k) 



1 

JTzy 
i 



Z 2 $ n (/c) + 2Z$ 12 (A;) + $ 22 (fc) 



$ u (k)-2$ 12 (k) + $ 22 (k) 



Z$ n {k) + (1 - Z)$ 12 (k) - $ 22 (fc) 



(19) 
(20) 



(1 + Z)" 

and the tildes mean a Fourier transform. 

A general study of the limit k — > of eq. (jlfij) would be tedious and we make the 
following simplifying assumptions: 

• The HS will all be of the same diameter cr Q = a. 

• We split the potential <& a p{k) into two parts, i.e. Q a p{k) = $^g(fc) + <& A {k) where 
^afl(^) > denotes the repulsive part and $ (A;) < denotes the attractive one. 
The attractive part Q A is necessary to induce a liquid- vapor transition in the RPA 
approximation. Since we have considered HS of equal diameters the repulsive 
potentials ^Jk) are used to mimic the soft core asymmetric repulsive interactions. 
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Moreover, at small k, we assume that 

*3»(*) = *5>(0) (l - {b% kf) + O(k') (21) 
^(jfe) = $ A (0) (l - (b A A;) 2 ) + 0{k A ) , (22) 

which implies that 

®AB(k) = $ AB (0) (1 - (b AB kf) + 0(A; 4 ), (23) 

where b A B are complicated functions of Z, 6^, 6 A , $^g(0) and <3> A (0) which will not be 
displayed here. The b% can be interpreted as the effective diameters for the ions. 
As a consequence of these assumptions we have a simple result 

G NN (k) G^ NQ (k) \ 

G NQ (k) G QQ (k) J 1 1 

( pr^r + ^NNik) P®N Q (k) ^ 

Gns{k) 

1 47T/?(7 2 

where Ghs(^) is the Fourier transform of the truncated two-body correlation function 
of a one-component HS fluid at the density p^. We are now in position to study eq. 
flUJ) in the limit jfe -»• 0. 

At = one finds that 

G NN (0) = - (25) 

1+/3^(0)Ghs(0) 

or 

GW(0) = - 1 , (26) 

z/!(p) +/3$ iV iv(0) 

where ^i(p) = <9z/ H s/<9p. 

It is worth noting that if $jV7v(0) > 0, then (jaw(0) is regular and there is no critical 
point (CP) in this case. If $^^(0) < 0, the isotherm compressibility \t = PGnn(0)/ P 2 
can diverge signaling the occurrence of a CP. It turns out ^3] that P\{p) is a positive 
convex function of the density with a single minimum at p c <r 3 = 0.248 with a value 
Vi, c = ^i(Pc) = 11.115. It follows from this remark that, in the RPA approximation, the 
critical density is p c and that the critical temperature is given by /3 c §nn(0) = —v\ iC . 
Hence, along the critical isochore and above T c , Gnn(0) behaves as 

G NN (0) = L (27) 

(P-Pc)*nn{0) 

yielding the critical exponent of compressibility 7 = 1 as expected. 

Following Aqua and Fisher jTSj we can cast the expressions for G^Nik) and Gqq(/c) 

as 

° NN{k) ~Wj ~^ikj~ (28) 

r (k) - S{k) 1 1 ~ S{k) (20) 
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m = i - ^Sgf- + o{k*) . (so 



- P&nnQWnn ~ HCK ^ JJ . (32) 



where Ai(fc) and A2(fc) are the eigenvalues of the matrix of the direct correlation functions 
Cab and 

WQ) fc 4 

lGir 2 q 

It is easily shown that at small k 

X 2 (k) = {(3- /3 C )$ NN (0) + ak 2 + 0(k A ) , (31) 

with 

1 d 2 c RS ( Pc ,k Q ~ (M Pc($nq(0)) 2 

a ~ 2 dk 

In the above eq., chs is the ordinary direct correlation function connected to Chs by 
Cns(k) = chs(^) _ 1/P- Therefore, along the critical isochore above T c we have 

G NN (k) = L — l -— for k -> , (33) 

(/3-/3 c )<W0)l + £ 2 ^ 

where the squared density- density correlation length is given by 

= ^ (34) 

(J3 - /3 e )$ NN {0) 

The positivity of £ 2 , i.e. the positivity of a, implies some restrictions on the various 
parameters of the model. As expected, the critical exponent of the correlation length 
is v = 1/2 and the Fisher exponent rj = 0. In the case where q = and $ fl = one 
recovers for a the result of ref. [Tl] . 

Finally the behavior of Xi(k) for A; — > is found to be 

\ 1 (k) = ^(l + ^-)+0(k*), (35) 



k 2 



K 



D 



where we have defined an effective squared Debye number as 

75* = £ , (36) 

1 + ^ QQ (0) P Z 

k 2 d = Anpj3q 2 Z being the usual squared Debye number. Hence, once again along the 
critical isochore we get for Gqq(A;) the expression, 

which is valid far and at the CP at small k. 

Some comments are in order. Firstly, away from the CP (/3 ^ f3 c ) we have 
Gqq(/c) ~ k 2 /4:Tc/3q 2 which means that the both SL rules are satisfied. This is in 
agreement with the results obtained in for the case when the direct correlation 
function is given by (jTBj). In the general case, at the CP (/3 = /3 C , p = p c ), = and 
we have for Gqq(/c) 

g «M-4^+%^W 4) - < 38) 
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As is seen from the above equation Gqq(A; = 0) = and thus the first SL rule is satisfied 
at the CP. By contrast, GQg(k ^ 0) is proportional to k 2 but with the coefficient that 
is not correct to ensure screening. Hence, the second SL rule is violated at the CP. 

However, in the special case of a symmetric model such that $f Q = $^jg and thus 
^NQ(k) e0 we have the simple result 

GQQW = TT^-T ( for allk ) > ( 39 ) 

and, in this case, the second SL rule is satisfied even at the CP. 

Our last comment will be for the charge- density correlation function, which along 
the critical isochore above T c behaves as 

s ^ = -WwTF*) +oihl) - (40) 

at small k. Therefore, for the asymmetric model G^Qik) ~ k 2 away from the CP when 
k — > and G NQ {k = 0) = -^ NQ (0)/(ATTq 2 a) ^ at the CP. Note that for the symmetric 
model Gnq = 0. 

In summary, we have studied density-density, charge-charge and charge-density 
correlation functions for the charge- and size-asymmetric model far and at the gas- 
liquid CP within the framework of the RPA and have shown that in a general case the 
second-moment SL rule is satisfied away from CP and is not satisfied at the CP. In the 
particular case of the symmetrical model the both SL rules are satisfied at the CP. It is 
worth noting that a more comprehensive study of the criticality in ionic fluids requires 
the fluctuations to be taken into consideration at the level higher than the RPA. This 
issue as well as a detailed analysis of an effect of the charge- and size- asymmetry on the 
phase behaviour will be presented elsewhere. 

References 



[l] 
[2] 

m 

[4] 
[5] 
[6] 
[7] 
[8] 
[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 
[17 



Weingartner H and Schroer W 2001 Adv. Chem. Phys. 116 1 
Fisher M E 1994 J. Stat. Phys. 75 1 
Stell G 1995 J. Stat. Phys. 78 197 

Stell G, Wu K C and Larsen B 1976 Phys. Rev. Lett. 37 1369 

Yan Q and de Pablo J J 1999 J.Chem.Phys. Ill 9509 

Panagiotopoulos A Z 2002 J. Chem. Phys. 116 3007 

Caillol J-M, Levesque D and Weis J- J 2002 J. Chem. Phys. 116 10794 

Luijten E, Fisher M E and Panagiotopoulos A Z 2002 Phys. Rev. Lett. 88 185701 

Patsahan O V and Mryglod I M 2004 J. Phys.: Condens. Matter 16 L235 

Netz R R and Orland H 1999 Europhys. Lett. 45 726 

Caillol J-M 2005 Mol. Phys. 103 1271 

Aqua J-N, Banerjee S and Fisher M E, arXiv:cond-mat/0410692 



Aqua J-N and Fisher M E 2004 Phys. Rev. Lett. 92 135702 
Caillol J-M 2003 Mol. Phys. 101 1617 
Yan Q and de Pablo J J 2002 J.Chem.Phys. 116 2967 
Panagiotopoulos A Z and Fisher M E 2002 Phys. Rev. Lett. 88 045701 
Evans R and Sluckin T J 1980 Mol. Phys. 40 413 



